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The influence of local order on the disordering scenario of small Wigner islands is discussed. A 
first disordering step is put in evidence by the time correlation functions and is linked to individual 
excitations resulting in configuration transitions, which are very sensitive to the local symmetries. 
This is followed by two other transitions, corresponding to orthoradial and radial diffusion, for which 
both individual and collective excitations play a significant role. Finally, we show that, contrary 
to large systems, the focus that is commonly made on collective excitations for such small systems 
through the Lindemann criterion has to be made carefully in order to clearly identify the relative 
contributions in the whole disordering process. 
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Many efforts have been intended in past years to un- 
derstand properties of mesoscopic devices in which inter- 
acting particles are confined. For instance, these parti- 
cles can be vortices in mesoscopic shaped superconduc- 
tors [l| , electrons in quantum dots [|[ , strongly coupled 
rf dusty plasma [3, trapped cooled ions vortices in 
superfluid HeA [6|, electron dimples on a liquid helium 
surface Q, vortices in a Bose-Einstein condensate [1] or 
colloidal particles @. More recently, in order to take 
advantage of the macroscopic scale to explore the prop- 
erties of such systems, we have proposed a new macro- 
scopic system consisting of N interacting charged balls 
of millimetric size free to move on a plane conductor and 
confined electrostatically, the temperature being simu- 
lated by a mechanical shaking [l(| ■ Using this system we 
observed the equilibrium configurations of the Wigner is- 
lands obtained for circular II il l and elliptic confinements 

m. 

Our previous study, essentially focused on static prop- 
erties, conclude that at low temperature and for a circular 
confining potential, the observed small islands (N < 40) 
present self-organized patterns constituted by concentric 
shells in which the balls are located. As it was widely 
discussed in the literature 13] [14] 15], this peculiar struc- 
ture is due to the competition between the ordering into 
a triangular lattice symmetry, which appears for infinite 
two-dimensional electrostatic systems, and the circular 
symmetry imposed by the confining potential. In the 
following, the Wigner islands will be described by the 
configuration (Nq — Ni — N2 — ..) where N% is the number 
of balls in the i th shell from the center. 

Surprisingly, the influence of the thermal fluctuations 
on the phase behavior which was an important question 
largely studied for the two-dimensional extended systems 
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has been little discussed in such small systems. The stud- 
ies devoted to this question are essentially numerical and 
focused on the temporal stabilities of the different config- 
urations [l6| [13], their spectral properties l|| 19[ or the 
mean displacements of the particles [l4j]. Let us how- 
ever indicate an interesting experimental work concern- 
ing the "mel ting " p rocess of colloidal interacting parti- 
cles islands [2fJ [21J . "Melting" in two-dimensional quan- 
tum electron clusters has also been recently studied in 

Phase transitions in two-dimensional large crystals 
have mainly been described by the changes in the asymp- 
totic behavior of spatial correlation functions. For in- 
stance, the Kosterlitz-Thouless-Halperin-Nelson- Young 
(KTHNY) theory predicts a two-step melting scenario ac- 
cording to which the liquid phase is reached when bond- 
orientational correlations become short-range [23| [2^] . 
Parallely, it has been suggested that temporal correla- 
tions should have the same behavior as the spatial ones 
[25I ] . This has recently been put in evidence experimen- 
tally for colloidal systems @ [26| . 

On the other hand, the previous studies on small in- 
teracting systems always refer to a generalization of the 
well-known Lindemann's method employed to describe 
the order-disorder transition for large systems [2?| , which 
considers the liquid phase is reached when the mean 
square displacements relatively to the lattice parameter 
of the particles go beyond the value 7a/ = 0.1 (which 
seems to be independent of the interaction [27| [HI), or 
0.05 for each coordinate. Note that it has been shown 
that, for infinite systems, the transition temperature ex- 
hibited there is the same as the one given by the correla- 
tion functions 29]. Let us underline that no such result 
is known for small systems. 

For the latter, the studies lying on Lindemann's crite- 
rion predict or observe that the shell-structured islands 
become less ordered while the temperature increases. 
They describe a two-step process corresponding to two 
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different transition temperatures. At very iow temper- 
ature, each particle is thermally excited in its local po- 
tential; a first transition appears at the temperature To 
when the orientational order between the shells is lost. 
This first transition is followed by a second one, at the 
temperature Tr, which corresponds to the emergence of 
the radial diffusion of particles between the shells, as well 
as an angular diffusion in the shells. For higher tem- 
perature, the initial order is completely destroyed. The 
transition temperatures To and are respectively iden- 
tified as the temperatures at which the intershell angular 
and radial mean square displacements have a rapid and 
strong increase [T3 |. 

In this paper, we discuss the influence of the local or- 
der on the "melting" . Indeed, for small islands, singular 
events such as a unique jump of only one particle from 
one shell to another one result to an important modifi- 
cation of the system since such an intershell jump corre- 
sponds to a transition between well-separated stable and 
metastable equilibrium configurations of different geom- 
etry. Consequently, such jumps will modify the bond- 
orientational correlation functions. We shall show that, 
for a given temperature, those configurations switches 
and thus the disordering process are controlled by the 
local geometry. In order to distinguish between those 
events and the collective excitations of the particles, we 
calculate the Lindemann parameters on each shell and 
for each configuration separately, without taking into ac- 
count jumps between shells. We determine the tempera- 
tures at which the collective excitations appear and de- 
stroy the ordered configurations. This unusual Linde- 
mann procedure is in agreement with our definition of an 
ordered system as a system always close to its most sym- 
metrical configuration. For instance, an intershell jump 
which does not destroy the local symmetry does not have 
to be considered as relevant in our description. This will 
explain why the obtained critical temperatures will be 
higher than the one determined by the usual Lindemann 
method which includes all kind of displacements in the 
calculation of the parameters, as in [14j . Notice that this 
question is specific to small systems since in large systems 
there is a continuum of states and individual excitations 
are hidden by collective ones. 

In order to evaluate the relative importance of these 
different contributions to the "melting" , we have exper- 
imentally observed the evolution of macroscopic Wigner 
crystals while the effective temperature is increased. To 
emphasize the contribution of the individual excitations 
with respect to the collective ones, we have selected a 
set of systems which have very different local symmetry 
for a similar number of balls in order to present different 
configuration transition behaviors for almost the same 
kind of collective excitations. We chose Wigner islands 
consisting in N — 18, 19 and 20 interacting particles con- 
fined in a circular frame. These systems, in spite of a 
very close number of balls, are very different from the 
local symmetry point of view and the resulting excita- 
tion energy spectra. Indeed, the "magic number" system 



N = 19 exhibits a three fold symmetry, as its ground 
configuration is (1-6-12). In fact, the latter shell can be 
divided into two subshells of 6 balls, as it was numeri- 
cally shown in [l3|, however, the difference between the 
two radii being rapidly of the same order as the thermally 
induced radial fluctuations, we will still refer to this con- 
figuration as the (1-6-12) one. By contrast, the ground 
configurations of N = 18 and N = 20 systems are consti- 
tuted by incommensurable shells (in the sense that the 
ratio between the number of balls on each shell is not an 
integer), their ground states being respectively (1-6-11) 
and (1-6-13). When the temperature increases, the two 
lower excited states of each system can be reached. In 
spite of the radial displacements due to the temperature, 
the ring-like structure remains and the configurations can 
be very well identified : 

• For 18 balls : (1-6-11), (1-5-12), (0-6-12) 

• For 19 balls : (1-6-12), (1-7-11), (1-5-13) 

• For 20 balls : (1-6-13), (1-7-12), (2-6-12) 

We can notice that only one ball jump from a shell to 
another is necessary to induce a configuration switch. 

These configurations correspond exactly to those com- 
puted in (l3j for logarithmic interparticle interaction po- 
tential, strongly suggesting this kind of interaction be- 
tween the balls, at least within the range of our exper- 
imental inter-particle distance [llj . This conclusion has 
been confirmed later by comparing the ground configura- 
tion obtained for elliptic confinement with the configura- 
tion of vortices calculated in similar shaped mesoscopic 
superconductors for which the inter-vortices interaction 
is logarithmic in the considered range of inter-particle dis- 
tance 12 1 [3(|. From the energetic point of view, the local 
symmetry differences between the various systems induce 
strong differences in their excitation energy spectra [l3j |. 
The configuration energies for N = 19 and N = 18 are 
well separated, however the gap between the ground and 
the first excited state is larger for the magic number sys- 
tem. By contrast, the latter is very small for N = 20. 

In section HI we present the experiments which vali- 
date the mechanical shaking as an effective thermody- 
namic temperature and we describe the parameters used 
to characterize the disorder, the configuration transitions 
and the collective excitations. The evolutions with tem- 
perature of these parameters for N = 18, 19 and 20 is- 
lands will be described and discussed in details in section 
ITIl In section IIII1 the respective influence on the disor- 
dering of the individual and collective excitations will be 
discussed. We will show that, according to the consid- 
ered parameter, different transition temperatures can be 
identified. In particular, the systems present an impor- 
tant configuration transition activity inducing disorder 
at a temperature smaller than those characterized by the 
Lindemann criterion. The "melting" will have to be con- 
sidered rather as a disordering than a real melting. 
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I. EXPERIMENTAL PROCEDURES AND 
CHARACTERIZATION PARAMETERS 

Our Wigner islands are constituted by millimetric 
stainless steel balls (of diameter d = 0.8 mm and weight 
m = 2.15 mg) located on the bottom electrode of a hor- 
izontal plane capacitor (a doped silicon wafer whereas 
the top electrode is a transparent conducting glass). An 
isolated metallic circular frame of diameter D = 10 mm 
and height h = 1.5 mm intercalated between the two elec- 
trodes confines the balls [11]. When a potential V is ap- 
plied to the top electrode (the bottom one and the frame 
being linked to the ground), the balls become monodis- 
persely charged, repel each other and spread throughout 
the whole available space. For the currently used po- 
tential V (around a few hundred of volts), the charge of 
each ball has been evaluated to about 10 9 electrons. The 
whole cell is fixed on a plate linked with two indepen- 
dent loudspeakers supplied by a white noise voltage. As 
we shall show later, this horizontal shaking results in an 
erratic movement of the balls which simulates an effective 
temperature that can be modified by tuning the shaking 
amplitude. In the experiment presented here, the system 
is always prepared in its ground configuration at low ef- 
fective temperature and the temperature is progressively 
increased. 

Throughout the experiment, images of the arrays of 
balls are recorded in real-time using a CCD camera, and 
their center of mass is detected. For each temperature, 
the positions of each ball are followed during 400 s, the 
minimum time between two snapshots being 100 ms. The 
characteristic time of the oscillation of a single ball in 
its local potential is, at low temperature, about half a 
second. Then, the choice of the total experiment time 
and the time interval allow us to track the balls (at least 
in the considered temperature range) and get statistically 
relevant data. The different configurations reached are 
then determined by counting the number of balls in each 
shell. In this procedure the radial limits of a shell are 
defined as the mean values between its radius and its 
neighbor's one, that have been measured in the ground 
state. They are independent from the temperature. 



A. Temperature calibration. 

Before studying the "melting" of such systems, strong 
attention has been paid to show that the cell shaking 
effectively results in a brownian motion of the balls al- 
lowing the identification of this shaking with an effec- 
tive temperature. We present here the experiments per- 
formed in order to calibrate this effective temperature 
and the "in situ" thermometer we have developed. 

As in [3l|, we used a system for which the energy is 
well-known : the calibration was obtained by the use of 
a single ball rolling on the silicon wafer which has been 
winded with an angle a = 25' from the horizontal plane. 
The ball can elastically bounce on a bottom wall. No 



electrostatic force is at stake and the only energy is the 
gravitational potential. For each voltage A applied on 
the loudspeakers, the horizontal distances x of the ball 
center from the bottom wall have been measured through 
the capture of a few thousand snapshots of the ball. The 
recorded random x positions are distributed following the 
density P[x). 

In order to compare this density with Bolzmann theory, 
it was fitted with the function 

, . mo tan a - E (*) 
P{x) = y e *bt 

where E{x) — mgx tan a is the potential energy of the 
ball and T the fit parameter which corresponds to the 
expected effective temperature, ks being the Boltzmann 
constant. 

Figure QJa) presents the experimental data obtained 
for the voltage amplitude A = 1.0 a.u. and the cor- 
responding fit. At evidence Boltzmann law is obeyed. 
This very good agreement being observed whatever the 
voltages A, we may conclude that the mechanical shak- 
ing corresponds to an effective temperature. This ef- 
fective temperature is obtained thanks to the fitting 
analysis, or more simply through the relationship T = 
mg < x > tana/fcs and a calibration curve T(A) has 
then been be obtained. As shown in figure QJb), the 
relation between A and T is affinc within the range re- 
quired to study the "melting" . Let us indicate that the 
temperature range is about 10 11 K, which has no other 
signification than the energy range. 

Finally, we have used this calibration to develop an"in- 
situ" thermometer which is constituted by a single ball 
trapped in a second circular frame located near the 
main one and submitted to the same voltage. The ef- 
fective temperature is determined by measuring the ra- 
dial mean square displacement < r 2 > of this unique 
ball for each given shaking amplitude A and by identi- 
fying this displacement with the temperature T(A) pre- 
viously determined by the calibration procedure. What- 
ever the various thermometer diameters tested in order 
to optimize the sensibility of this thermometer (D = 
5, 6, 7, 8 and 10 mm) and the applied potential V (V = 
700,800,900,1000 V) the mean square displacement < 
r 2 > varies linearly with the temperature. Same lin- 
earity is observed for the mean square speed. Figure [2] 
presents the variation < r 2 > (T) corresponding to the 
potential V — 900 V and the retained diameter D = 6 
mm. This thermometer will give a precise determination 
of the temperature for any experiment to come, indepen- 
dently from the variations associated to the total weight 
of the support or to the loudspeakers' ageing. 



B. Correlation functions 

The bond-orientational correlation function ge(t) is rel- 
evant in order to characterize phase transitions in two- 
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FIG. 1: Temperature calibration : (a) x position distribution for a trapped ball moving on an inclinated plane (inset) and 
submitted to mechanical shaking (log scale). The good fit with Boltzmann law validates the mechanical shaking as an effective 
temperature, (b) Calibration of the effective temperature. 
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FIG. 2: Mean square displacement < r 2 > of the ball plotted 
versus temperature for D=6 mm and V=900 V. 

dimensional systems. It is denned by 

g 6 (t) = | < e Me(t+t )-e(t )) > | 

where 9(t) is the angle of a fixed bond between two par- 
ticles and <> denotes an averaging over all bonds. In 
infinite two-dimensional systems, g§ tends to a constant 
roughly equal to 1 at low temperature. The system is 
then like an ordered crystal. If the temperature is higher 
than a temperature XJ, a strong decay with time is ob- 
served, which denotes a liquid phase. 

Such a dynamical criterion offers experimental facili- 
ties since it is often easier to have long time acquisition 
rather than to observe a large sample of particles. In par- 
ticular, we can readily use this criterion to characterize 
order-disorder transitions in small systems. 



Note that the factor 6 in the exponential is adapted 
for hexagonal lattices, therefore we can expect the limit 
value at low temperature to be lower than 1 for our small 
systems in which the three-fold symmetry is broken. On 
the other hand, the averaging being made only on N(N— 
l)/2 links, we can't expect ge to go until zero value in 
the liquid phase, even at large time. These points apart, 
this parameter still measures the correlations and we can 
expect to measure the same global behaviors according 
to the temperature. 

C. Configurations transition parameters 

Two parameters easy to get experimentally have been 
identified in order to characterize the configuration tran- 
sitions. The first of them is the jump rate Rj : 



where Nj(t) is the total number of configuration 
switches during the period t. This parameter is an in- 
dicator of the "transition activity" of the system. Quali- 
tatively, this parameter is small at low temperature when 
only a few transitions occurs and increases strongly with 
temperature when energetical barriers can be overcome. 

A second way to characterize more quantitatively the 
transition rate is to measure the mean time required to 
escape from each state or the mean residence time in 
each state. Let us consider for instance a Wigner island 
at low temperature ; this system can be understood as 
a two level system characterized by a ground state E g 
and one metastable state E m (the only one that is reach- 
able if the temperature is sufficiently low). The thermal 
fluctuations induce transitions between those two con- 
figurations. These transitions in the real space can be 
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mapped in the phase space by the jump of a Active par- 
ticle from a well to another, characterized by their depth 
Eg and E m and a saddle point E s . 

Under the assumption that all escape attempts are in- 
dependent and of weak probability, the probability for 
the particle to stay in a state of energy E during a time 
r is : 



P(t) = — e - T / T0 

TO 

where tq =< r > is the mean residence or Kramers 
time. It depends on the energy barrier and on the tem- 
perature and is given by [32| : 

< r >= T r e^- E ^' kBT 

where E s — E is the barrier energy, and r r is the relax- 
ation time within the well. 

The energetical barriers can be measured by the slope 
of the curve describing the variations of the mean resi- 
dence times in log scale as a function of 1/T. Similarly, 
the spectrum of the excitation energies is determined 
through the ratio of mean residence times in the two 
wells. 

Qualitatively, this analysis can be extended to the case 
of high temperature for which the system explores more 
than the two first levels and reaches higher excited lev- 
els. But the possibility for a particle to escape a well in 
different ways whose relative weight should depends on 
the geometry leads to a more acute problem. 

We have tested the validity of this analysis in the case 
of the two- level systems N = 5 and N — 6. These 
two cases are interesting since they involve in their re- 
spective configuration transitions the two kinds of in- 
dividual jump observed in the Wigner excitations. For 
N = 5, the ground configuration consists in a unique 
shell (5) whereas the metastable configuration (1-4) re- 
quires a centripetal displacement of a ball. By contrast, 
the N — 6 ground state is a centered configuration (1-5) 
and the metastable configuration (6) is reached after a 
centrifugal displacement of ball. 

In figure [3] we present the variation with the temper- 
ature of the ratio < T centered > / < T circ i e >, where 
Tcentered stands for the residence time in the centered 
configuration, respectively (1-4) or (1-5), and r C i,. c ; e de- 
notes the residence time in the one-shell configuration, 
respectively (5) or (6). For the two selected systems these 
variations obey Kramers' relation. 

According to this boltzmanian description, we will 
characterize in the following the configurations transi- 
tions by Rj and by the mean residence times in each 
configuration. These parameters depend on the config- 
uration energy spectra and thus, are extremely sensitive 
to the local symmetry of the configurations. 
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FIG. 3: Evolution with the temperature of the residence mean 
times ratio for two-level systems : 5 balls (■) and 6 balls (•). 
Ratios are in log scale. Note the signs of the slopes that 
prove that (5) and (1-5) are the ground states for the 5-ball 
and 6-ball systems respectively. 

D. Lindemann-like criterion 

In order to explore the collective excitations of Wigner 
islands through a Lindcmann criterion, the mean square 
deviations of the balls from their equilibrium locations 
have to be calculated. 

For three-shell configurations and circular symmetry 
[Til ], we have to calculate, for each configuration and for 
shells 1 and 2, the radial displacements 

s i— 1 

the relative angular intrashell displacements 
1 Ns 

< = W^ {< {dl - > - < °* °n > 2 )/^o, 

s i—l 

and, for each configuration and for shell 1, the relative 
angular intershell displacements 

1 Ns 

"m = £(< V* - B **¥ >-<o<- K > 2 )/o 2 

i—l 

where are the polar coordinates of a ball with 

respect to the center of the confining frame, i\ indicates 
the right neighbor of the ball i, and 12 indicates its nearest 
neighbor in the surrounding shell, which is determined 
every snapshot, ro = 1 /s/nn, where n is the balls density, 
is the mean radial free space for a ball and 9q — 2ir /N s is 
the mean interball angular distance in the shell consisting 
in N s balls. <> indicates an average over time. 

These deviations are relative displacements (with re- 
spect to another ball or the center of the system), there- 
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fore they are relevant in order to exhibit a Lindemann- 
like criterion comparable to the one used for larger two- 
dimensional systems [13] ■ However, we will not be able 
to discuss as in the latter article about the values of the 
dimensionless parameter T = E/ksT, where E is a typi- 
cal interaction energy between two particles, since we do 
not have numerical values for the potential energies in 
our system at the present time. 

Note that we have discriminated between the differ- 
ent configurations. This is allowed since the radial dis- 
placements, as we shall show, are around a sixth of the 
distance between the two shells whatever the effective 
temperature, then the shells are always well identified. 
This choice is motivated by our requirements of precise 
information about the influence of the local order on the 
"melting" . The temperature dependencies of the dis- 
placements averaged over all the configurations have also 
been calculated. This procedure is not the same as the 
one used in where the distinction between shells is 
only made at the beginning of the numerical simulation, 
so we shall expect lower values for our radial displace- 
ments that do not include intershell jumps. 



II. IDENTIFICATION OF TRANSITION 
TEMPERATURES 

A. Bond-orientational decorrelations 

On figure HJa-c) we present typical time correlation 
functions ge(t) for five different temperatures and for the 
three systems, for which the same behaviors are observed 
: at low temperature, a constant value is reached, indi- 
cating an ordered state. At higher temperature, the cor- 
relations decay strongly. As in large systems, we can de- 
fine a "melting" temperature Ti (N). We notice that the 
N = 19 system keeps an ordered structure denoted by a 
quasi-constant correlation function until a much higher 
temperature than the two other systems. Indeed, the 
transition temperatures T/(iV) are sensitively different 
according to the number of balls N : T/(18) is less than 
13 x 10 n K, whereas T/(19) is more than 18 x lO n K and 
T/(20) between both. Note our goal is not to measure 
precise transition temperatures, but to put in evidence 
the mechanisms involved in the disordering. As we shall 
see, the different transition temperatures are sufficiently 
separated to do such an analysis. In the following, we 
will then investigate the two expected mechanisms for 
this disordering : configuration transitions and collective 
excitations through Lindemann criterion. 



B. Configurations transition 

The first indication of the configuration transition ac- 
tivity can be evidenced by observing the evolution of the 
jump rate Rj while the effective temperature increases. 



These variations are shown on figure[4jd) for the different 
studied systems. 

Whatever the number of balls, Rj presents the same 
qualitative behavior: it is close to zero at low tempera- 
ture and increases strongly at higher temperature. These 
variations correspond to a progressive augmentation of 
the number of transitions activated at the effective tem- 
perature. Whatever the temperature, the smallest jump 
rate is associated to N — 19 system and the one asso- 
ciated to N — 20 is always smaller than the one corre- 
sponding to N = 18. In order to describe more quanti- 
tatively these behaviors and their differences, let us in- 
troduce transition temperatures which characterize the 
"beginning" of the Rj increases. We have chosen to 
name transition temperature the temperature at which 
the ground configuration begin to switch. Let us nev- 
ertheless indicate that this transition temperatures does 
not correspond to an actual sudden transition since the 
jump rate Rj rises progressively. Note that experimen- 
tally, the infinite time limit in the Rj definition corre- 
sponds to the maximum measurement time and its finite 
value could alter the Rj value. However a complete anal- 
ysis has shown that the variation with temperature of Rj 
is independent from this experimental limit provided that 
this limit was largely higher than the mean residence time 
of the system in each configuration. In our experiments 
we chose t max = 400s which satisfied this condition. 

The analysis of the evolution of the residence times dis- 
tribution shows that, for N = 19, the transition tempera- 
ture Tj(19) is 20 x 10 11 K, much higher than the tempera- 
tures Tj(18) and Tj(20) which are respectively 14x 10 n K 
and 17x 10 11 K. At larger temperature, the second excited 
states are reached. The corresponding temperatures are 
respectively 21 x 10 n K, 16 x 10 n K and 21 x 10 n K for 
N = 19, 18 and 20 and are indicated by arrows on figure 
I2£d). Those temperatures might correspond to the crit- 
ical temperatures given in [161 ] where the rate of radial 
jumps is calculated. This is exactly the same parameter 
as our jump rate since each configuration switch involves 
only one ball jump from one shell to another. 

These behavior differences between the systems are 
still more obvious when we study their mean residence 
times in the ground state. On figure [5] we have plotted 
the logarithm of these times < t > versus the inverse of 
the effective temperature. We obtain a linear variation, 
which is in agreement with the Boltzmann law descrip- 
tion presented above, the curves slopes being equal to 
the barrier heights for escaping the configuration. Wc 
can observe that these slopes are different according to 
the number of balls. The highest slope is associated to 
the ground state N = 19 and is equal to 192 x 10 11 K, 
those corresponding to N — 18 and 20 being respectively 
equal to 99 x 10 n K and 107 x 10 n K. This indicates that 
the ground state for N = 19 is much more stable than 
the ground states associated to TV = 18 and 20, whose 
barrier heights are very similar since their corresponding 
slopes are roughly identical. 

This analysis can be completed by studying the excited 
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FIG. 4: Variations of orientational correlation function and jump rate with temperature. Bond-orientational time correlation 
function for five temperatures : (■) T = 9.6 x 10 n A", (•) T = 13.0 x lO 11 ^, (A) T = 17.8 x lO 11 ^, (T) T = 23.0 x 10 n i\, 
(♦) T = 30.2 x 10 n i\ ; (a) N = 18, (b) N = 19 and (c) N = 20. (d) Jump rate versus temperature. Arrows show the 
temperature at which the second excited state is reached. 



states (figureE). The slopes for iV = 18 is 64 x 10 n K and 
is equal to 101 x 10 11 K for N = 19 ; in both systems, this 
suggests that the barrier is much lower than the ground 
state's one. On the contrary, the slope associated to the 
first excited state of N — 20 is equal to 96 x 10 11 K, al- 
most the same as the barrier for the ground state. All 
these results are ingood agreement with the energy lev- 
els calculated in [131 ]. Let us indicate that even when 
they are reachable, the other excited states are not easy 
to study since their higher energy involves too small sta- 
tistical occurrences. Finally, these measurements allow 
us to determine the first excitation energy associated to 
each system. These energies are 35 x 10 11 K, 91 x 10 11 K, 
11 x 10 n K respectively for N = 18, 19, 20. 

These measurements show that the "magic number 
" system corresponds to the deeper ground state and con- 
firms its strong stability in comparison with the two other 
systems. As expected, the commensurability influences 
strongly the depth of the well, and consequently the tran- 
sition temperature in Rj. 

Let us conclude this section by indicating that obvi- 
ously the jump rate Rj is related to the mean residence 
times. For instance, for a two-level system, Rj is simply 
equal to 2/(< ri > + < >), where the subscripts 1 



and 2 stand for the two levels. This relation is satisfied 
and confirms the self consistency of our results, at least 
up to a temperature at which a third level is reached. 



C. Mean square displacements 



We now turn to the study of the "melting" through 
Lindemann-like criterion. We will successively present 
the radial, intrashell and intershcll mean square displace- 
ments averaged over all the configurations. In order to 
explore more precisely the relation between the local or- 
der and these collective parameters, we have compared 
them to the corresponding parameters before averaging. 
We shall show in particular that the procedure of averag- 
ing over the configurations, commonly used in literature, 
mask actually subtle effects resulting from the configura- 
tions transitions, even though they don't infer that much 
on transition temperatures. 
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FIG. 5: Evolution with the temperature of the residence mean 
times in the ground state (the 400s limit corresponds to the 
recording time but it is actually infinite). Times are in log 
scale. 
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FIG. 6: Evolution with the temperature of the residence mean 
times in the first excited state. Times are in log scale. 



1. Radial displacements 

The temperature dependencies of the radial mean dis- 
placements averaged over all the configurations are pre- 
sented for the three different numbers of balls on figure 
Ufa). 

The displacements corresponding to the inner shell 
(shell 1) present similar behaviors for the three systems. 
They vary regularly from 0.005 at low temperature and 
go on increasing until 0.05 at the highest experimen- 
tal temperature, the highest and the lowest radial dis- 
placements being respectively associated to the systems 
N = 20 and N = 19, whatever the temperature. In 
this experimental temperature range, the low values of 
u^, smaller than the Lindemann criterion, indicate that 
each shell remains very well identified and that we can 
always discuss qualitatively the results in terms of shell 



entities. Notice that the highest obtained temperature 
corresponds to the ball-tracking limit, but we could fore- 
cast that the value would increase strongly beyond 0.05 
after this temperature limit, since this value is far from 
being the maximal possible value, even though jumps 
from one shell to another are not taken into account. 
From this point of view, we can consider the highest ex- 
perimental temperature is very close to the radial tran- 
sition temperature transition Tr of these systems. The 
radial displacements of the outer shell (shell 2) are quali- 
tatively similar to those observed for the shell 1 although 
their values are smaller and vary from 0.005 until 0.03. 
This can be explained by the fact that this shell is sub- 
mitted to a regular and constant potential resulting from 
the confinement frame whereas shell 1 is submitted to 
fluctuant potential from both sides [lij • 

We have seen that geometrical considerations are re- 
sponsible for the different configuration transition rates 
observed. It is therefore natural to examine if it is the 
same for the displacements. The radial mean displace- 
ments for the different configurations of the different sys- 
tems are presented on figure [5J Notice that data for the 
excited states with low residence time present higher sta- 
tistical error due to the smaller number of their occur- 
rences; this is the case for instance for the second excited 
state for 18 and 20 balls and for both excited states for 
19 balls. 

As for the averaged curves, the radial displacements 
in the two shells of the systems in their ground con- 
figurations present regular increases with temperature. 
More precisely, for shell 1, we can observe that the 
N = 19 (1/6/12) and N = 20 (1/6/13) curves are iden- 
tical whereas the displacement corresponding to N = 18 
(1-6-11) is higher whatever the temperature. For shell 2 
the displacements are identical for N = 18 (1/6/11) and 
N = 19 (1/6/12) whereas those associated to N = 20 
(1-6-13) presents a higher value. The temperature de- 
pendencies observed for the excited states look like those 
associated to the ground states, with a constant switch. 
The highest v% values are now observed for TV = 20 (1-7- 
12) in the case of the shell 1 and for N = 18 (1-5-12) in 
the case of the shell 2. 

We cannot define a precise rule to explain the different 
behaviors. It seems however that, for a given tempera- 
ture, the more balls in the shell, the larger the radial dis- 
placements. In the case of equality (that's to say, for the 
first shell, the couples (1-6-11), (0-6-12) and (1-6-13), (2-6- 
12) and for the second shell the couples (1-5-12), (0-6-12) 
and (1-7-12), (2-6-12)), we observe that the commensu- 
rable states have lower radial displacements. Those dif- 
ferences are higher than a possible distortion due to a 
renormalization by a different N s . 

The first statement can be simply understood by the 
fact that orthoradialy squeezed balls have more kinetic 
energy to involve in their radial movement. In the case 
of a commensurable state, all the balls (for the first shell) 
or half of them (for the second shell) find themselves in 
front of a repelling ball of the other shell, which restrains 
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FIG. 7: Mean square displacements averaged over the three configurations versus temperature : (a) radial displacements v% 
for both shells, (b) intershell angular displacements Ug 2 , (c) intrashell angular displacements u|j in shell 1 and (d) intrashell 
angular displacements Ugi in shell 2 



their radial fluctuations. Note this result is opposite to 
what was reported for quantum dots in [22] : the au- 
thors exhibit higher radial displacements for magic num- 
ber N — 19. 

Beyond these specific results, the comparison between 
the averaged and non averaged radial displacements show 
that even if some differences can be exhibited concerning 
their precise and relative values, global behavior remains 
the same : each shell will radialy melt at the same tem- 
perature whatever the configuration. 



slowly, the angular ones begin to vary linearly and change 
very rapidly at almost the same temperature whatever 
the number of balls. By contrast, this similarity of behav- 
iors is not observed for the outer shell : whereas the intra- 
shell orthoradial displacements are identical for N = 18 
and N = 20, in the case TV = 19, it remains smaller than 
0.005 and without rapid rise. 



2. Intrashell angular displacements 

In the case of the angular mean displacements, the av- 
eraging does not introduce any distortion in the analysis, 
as it can be seen on figure [5] So we will only discuss 
the mean displacements averaged over the configurations 
(figures [7tc) and (d)). Like the radial displacements, 
the intrashell orthoradial displacements associated with 
the inner shell of ground configurations increase with the 
temperature. However, whereas the former keep on rising 



The changes are observed when the displacements 
reach the critical value 0.05 in accordance with Lindc- 
mann criterion. Thus, we can define from these data a 
transition temperature To or, rather, a temperature in- 
terval centered temperature on To (between 28 x 10 11 K 
and 30 x 10 11 K) after which intra shell oricntational or- 
der is lost in a given configuration for N = 18 and 
N = 20. We can expect that the corresponding tem- 
perature is not far from the experimental limit in the 
case of TV = 19 since the balls in the shell 1 have begun 
to be non-correlated whereas the shell 2 is still rigid. 
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3. Intershell angular displacements 

Let us now consider the intershell relative displace- 
ments which measures the ability of the two shells to 
find a stable position one with respect to the other. Their 
variations with temperature averaged over configurations 
are presented in figure[7jb) for the three systems. We can 
notice that even at low temperature these displacements 
are high and until T w 18 x 10 11 K, the intershell angu- 
lar displacements are larger than the corresponding in- 
trashell and radial displacements. This indicates that at 
low temperature these intershell movements are the main 
effects resulting from the thermal fluctuations. Moreover, 
this temperature being smaller than To and T^, we are 
allowed to discuss these results in terms of well defined 
rigid shells in which the balls are regularly spaced. 

At large temperature, T > 18 x 10 11 K, the inter-shell 
angular displacements reach the same finite value what- 
ever the system, this value being in fact its theoretical 
maximum. In this temperature range, the thermal en- 
ergy is sufficient to overcome the barrier energies, which 
correlate the inner and outer shells. Since disorder is 
characterized by a deviation of the whole island from the 
symmetrical situation, note that the intershell displace- 
ments are calculated considering the nearest neighbor at 
each time step, in order to take into account the invari- 
ance towards some rotations. Consequently, we will al- 
ways have a small maximum value and will never observe 
a steep rise as in flij . 

Below this temperature, a temperature range in which 
the effect of the local order has an essential play is well 
exhibited. In this temperature domain, the systems are 
mainly in their ground state and the intershell angular 
displacement is lower for 19 balls than for 18 and 20 
balls, the latter remaining of the same order. This can 
be simply explained by commensurability arguments : 
let us consider the shells as rigid rings ; in the 1-6-12 
configuration, shell 1 is submitted to a 27r/12 periodic 
potential due to shell 2, then each of its balls can find 
itself in a potential well. On the contrary, in any other 
configuration uncommensurability implies that, if rigid, 
the two shells can't find a position for which every ball 
will be located in a minimum of energy, hence a higher 
instability. 

III. DISORDERING SCENARIO 

The study of the transition temperatures Tj through 
correlation functions show that they depend strongly on 
the local order. Indeed, those temperatures are very dif- 
ferent for the three systems, the N — 19 "magic system" 
being more stable. As for temperatures around Tj , the 
systems can still be seen as sets of shells, comparison can 
be made between the temperatures Tj and the tempera- 
tures Tj of first configuration transition. Since they are 
very similar, we can infer that configuration transitions 
play an important role in the disordering, at least for in- 



termediate temperatures. It can be simply understood 
by the fact that the islands are invariant by many rota- 
tions whatever the configuration, then a cycle of configu- 
ration transitions can induce disorder (from the correla- 
tion point of view), since the initial particles positions are 
not necessary recovered at the end of the transition cy- 
cle. Moreover a configuration switch implies a complete 
reorganization in the shells. This disorder mechanism in- 
volving only one particle jump cannot be evidenced by 
the global parameters as the mean displacements used in 
the Lindcmann approach. Note that even at low temper- 
ature, the orientational order between adjacent shells is 
already lost while their internal order is conserved ; this 
suggests that the trajectory of the jumping particle might 
be not only radial but also orthoradial, taking advantage 
of the relative rotation of the two implied shells. This 
is also in accordance with the fact that commensurabil- 
ity considerations play a role for the intershell rotations 
as well as for the height of the energetical barriers. On 
the other hand, mere intershell rotation is not sufficient 
to induce disorder as defined by the correlation function, 
since the system will periodically find itself in the same 
position as the initial one. 

If we focus now on collective displacements, a first tran- 
sition occurs at the temperature To > Ti, which corre- 
sponds to the emergence of the angular intra-shell diffu- 
sion. This stage in which the systems can be considered 
as independent shells remains until a second transition 
temperature Tr involving diffusion of particles between 
the shells. At higher temperatures, the shell structure 
disappears. Figure [TU] shows the mean crystals (superpo- 
sition of the different positions of the particles) obtained 
for instance for N = 18 at three key temperatures. Fig- 
ure fTUTb) illustrates the role of the transition activity 
in the loss of the orientational order : according to the 
values of the angular displacements, particles positions 
should be distinguishable. 

This kind of two-stage "melting" had been suggested 
by previously cited works, especially in (l4| . where the 
authors focus on a global Lindcmann criterion that in- 
cludes different mecanisms that we have exhibited there. 
In particular, we have distinguished two contributions in 
the radial displacements, namely the individual jumps 
and the mean behavior. In colloidal systems 20] [21] 
the system is arranged at low temperature in a shell- 
like structure. It also exhibits a very similar behavior 
when the temperature increases, excepted a re-entrant 
ordered phase which was not observed in our case, this 
phase being specific to the hard wall confinement [33(. 

In [22J, the authors studied the melting of N = 19 
and N — 20 quantum particles interacting with coulom- 
bic interaction. They described the melting as a two 
stages process : first an orientational inter and intrashell 
disordering and then a radial melting at higher temper- 
atures. The relative positions of the transition tempera- 
tures found here are in good qualitative agreement with 
their results. In addition, we can clearly distinguish two 
phases in the angular disordering, namely an intershell 
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N Ti Tj To Tr 

18 T( < 13 14 29 T H > 30 

19 T > 18 20 30 T R > 30 

20 13 < T < 18 17 29 T H > 30 

TABLE I: Summary of the different transition temperatures 
(in lO 11 ^). 

rotation and then an intrashell melting. This last point 
is also presented in numerical works on classical coulom- 
bic particles [27| . 

Whereas local geometry have an influence on the con- 
figuration transitions through commensurability consid- 
erations, we have shown that its effects are neglectable for 
the intrashell displacements as well as for the radial ones. 
On the other hand, correlation functions define very dif- 
ferent temperatures for our three systems (all the tran- 
sition temperatures are summarized in table H} . More- 
over, and contrary to large systems, those temperatures 
are lower than the transition temperatures given by the 
Lindcmann criterion. We have then to consider that be- 
yond the well-known Lindemann scenario, there are other 
sources of disorder, namely the configuration switches. 

IV. CONCLUSION 

In this paper we show that for a system constituted 
of a small number of interacting particles like Wigner is- 
lands, an increase of temperature results in a disordering 
of the system more than a real melting. This disordering 
process is very sensitive to the local order of the explored 
configurations. 

This disordering results from both individual excita- 
tions that induce configuration transitions and collective 
excitations. This process is marked by three different 
transitions. The temporal correlation functions which 
describe the correlation loss of the system exhibits an 



exponential decrease at the "liquid transition tempera- 
ture" Ti (named after the usual convention). This first 
transition is identified without ambiguity as correspond- 
ing to the increase of configuration transitions between 
the stable and metastable states of the system. So the liq- 
uid transition temperature depends strongly on the local 
order as the transition rate. At larger temperatures, two 
other transitions appear, To and Tr characterizing re- 
spectively intrashell and intershell diffusion. These tran- 
sitions are evidenced by the change of the mean square 
displacements with temperature and correspond to col- 
lective excitations. The local order play a less significant 
role in the transition temperature values. 

This disordering process and the importance of the lo- 
cal order on these temperatures are due the small num- 
ber of configuration states explored by the system at a 
fixed temperature. From this point of view, it is a spe- 
cific characteristic of small systems. Indeed, these effects 
cannot be observed for large systems since the number of 
explored configurations is large enough to mask the in- 
dividual excitations in the collective ones, hence the co- 
incidence between the transition temperatures described 
by the correlation functions and those identified by the 
Lindemann criterion. 

Whereas the T transition was never discussed before, 
the two following successive transitions have previously 
been mentioned in the literature. The corresponding 
analyses are in agreement with our results but their ap- 
proach is very different. We show that the description of 
the transition from well organized arrays towards liquid 
state resulting from successive excitations requires more 
detailed analyses than the single use of Lindemann crite- 
rion. 

Finally, let us conclude by suggesting that small 
Wigner islands that we proposed could be good can- 
didates in order to explore experimentally the thermo- 
dynamic laws dedicated to small systems of interacting 
particles[34j|. 
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